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C/^ \ Abstract 

Cf^ ' We analyze the responses of a quantum search algorithm to an external monochro- 

. matic field and to the decoherences introduced through measurement processes. The 

\ external field in general affects the functioning of the search algorithm. However, 
depending on the values of the field parameters, there are zones where the algo- 

■ rithm continues to work with good efficiency. The effect of repeated measurements 

jy-^ , can be treated analytically and, in general, does not lead to drastic changes in the 

\ efficiency of the search algorithm. 
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. 1 Introduction 



In the real world the concept of "isolated system" is an abstraction and ideal- 
ization. It was constructed to help understand some phenomena displayed by 
real systems which may be regarded as approximately isolated. Since dissipa- 
tion is a macroscopic concept, there has been little interest in it during the 
initial development of Quantum Mechanics. But, since about 40 years dissipa- 
tion has been incorporated into the quantum description to make possible the 
understanding of processes such as ionization of atoms, radiation fields inside 
a cavity or simply the decoherence caused by the interaction between a system 

^ Corresponding author. E-mail address: auyuanet@fing.edu. uy 

^ Permanent address: Instituto de Fisica, Universidade Federal do Rio de Janeiro 

CP. 68528, 21941-972 Rio de Janeiro, Brazil 



Preprint submitted to Elsevier Preprint 



1 February 2008 



and its surroundings. The recent advances in technology that have made possi- 
ble to construct and preserve quantum states, have also opened the possibility 
of building quantum computing devices [1,2,3,4]. Therefore, the study of the 
dynamics of open quantum systems becomes relevant both for development of 
these technologies as well as for the algorithms that will run on those future 
quantum computers. 

Grover's quantum search algorithm [5], which locates a marked item in an 
unsorted list of N elements in a number of steps proportional to y/N, instead 
of the 0{N) of the classical counterpart, is one of the more studied. This search 
algorithm has also a continuous time version [6] that has been described as 
the analog analogue of the original Grover algorithm. 

We have recently developed a new way to generate a continuous time quantum 
search algorithm [7]. In that work we have built a search algorithm with con- 
tinuous time that finds a discrete eigenstate of a given Hamiltonian Hq, if its 
eigenenergy is given. This resonant algorithm behaves like Grovers's, and its 
efficiency depends on the spectral density of the Hamiltonian Hq. A connec- 
tion between the continuous and discrete time versions of the search algorithm 
was also established, and it was explicitly shown that such a quantum search 
algorithm is essentially a resonance between the initial and the searched state. 

Recently [8,9,10] the response of Grover's algorithm to decoherences has been 
analyzed. Here a similar study on our proposal for a quantum search algorithm 
is performed. We rapidly review this method in the following section, and 
then proceed to study how this resonant algorithm can be affected by some 
interactions with the environment, as follows. In section 3, we subject the 
system to a monochromatic external field, and, in section 4, decoherences are 
introduced by performing a series of measurements. Conclusions are drawn in 
the last section of this work. 



2 Resonant Algorithm 

Let us consider the normalized eigenstates {\n)} and eigenvalues {£„} for a 
Hamiltonian Hq. Consider a subset N of {\n)} formed by N states. Let us call 
\s) the unknown searched state in N whose energy Eg is given. We assume it is 
the only state in N with that value of the energy, so knowing Eg is equivalent 
to ''marking" the searched state in Grover's algorithm. 

In the resonant quantum search algorithm a potential V, that produces the 
coupling between the initial and the searched states, is defined as [7] 



y{t) = \P) 01 exp {iusjt) + \j) {p\ exp {-iujsjt) , 
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where the eigenstate \j), with eigenvalue ej, is the initial state of the system. 
This initial state is chosen so as not to belong to the subset N. Above, \p) = 
^ |n) is an unitary vector which can be interpreted as the average of the 

set of vectors in N, and cUsj = Sj — Sg- 

The objective of the algorithm is to find the eigenvector |s) whose transition 
energy from the initial state \j) is the Bohr frequency ujsj- In order to perform 
this task, the Schrodinger's equation, with the Hamiltonian H = Hq + V{t), 
is solved. The wavefunction, is expressed as an expansion in the eigen- 

states {\n)} of Hq, |^(^)) = Z^m '^m(^) ^xp (— i£j„t) \m). The time dependent 
coefficients amit) have initial conditions aj(0) = 1, am{^) = for all m ^ j. 
We take units such that h — 1. 

After solving this equation, the probability distribution results in, 

Pj ~ cos^in t) , 

Ps ~ sm^{n t) , (2) 
Pn — 0, n ^ j and n ^ s, 

where fl = Prom these equations it is clear that for r = ^ a measurement 
has a probability of yielding the searched state very close to one. This approach 
is valid as long as all the Bohr frequencies satisfy cUnm ^ ^ and, in this case, 
our method behaves qualitatively like Grover's. 



3 Interaction with an external field 



The previous results have shown that the resonant search algorithm can be re- 
duced to essentially a quantum system with two energy levels. This reminds us 
of the quantum optics problem of a two-level atom driven by a radiation field, 
and the question arises on how the algorithm behaves when an external field 
is interacting with it. Let us then consider the time-dependent Schrodinger 
equation 

i^^-[Ho + v{t) + r{t)]mt)), (3) 

where the external monochromatic driving potential field 

r(i) = rosina;o^, (4) 

where Fq and ujq are the field's amplitude and frequency, respectively. Re- 
placing the expansion for in the equation above, we obtain the set of 
equations for the amplitudes am{t) 

^^^^-iY,{n\V{t) + r{t)\m)am{t)e^p{-icUnmt) , (5) 

at -tr 
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where ujrun = " ^16 Bohi frequencies. Inserting the definitions (1) and 
(4) into (5), we find, 

^""f^^ = -i sin uot V r„„a^(i) exp {-iunmt) , (6) 
at tr 



for n ^ N and n ^ j, and 

danjt) i_i 

dt ~ y/N\ 



[1 - 6nj) aj{t) exp [i {ujn + uisj) A 

+Snj ^rnit) exp [-i {ujjm + (^sj) t] } (7) 



a n or n — j. Above = (n|ro|m). 

As in the general time- dependent perturbation problem, we are in the presence 
of several time scales involved in the process. One is the fast scale associated to 
the Bohr frequencies, uJnm, another is the slow scale connected to the frequency 
Q — 1/ y/N, related to the natural oscilations of the system in the absence of an 
external field. There is also a third time scale, associated to the frequency ujq 
of the external field T{t), which in principle could take any value. We consider 
here cases where ujq < Q. In this case the presence of r(t) is relevant, so we 
do not make a restriction on the range of validity of the algorithm. 

Following the general procedure of time-dependent perturbation theory, we 
integrate the previous equations over a time interval much greater than the one 
associated to the fast scale. In this way, terms having small phases dominate 
and the others average to zero. Therefore the equations become 



1 

-Msit) — iTjjaj{t) sina;ot , 



daj{t) 

~dr-~7N'^'^'''~'^"^'' 

^^jf^ ~ l=aj{t) - iTssttsit) sinujot , (8) 

clt Y W 

- -i^nndnit) sina;ot for all s,j . 

These equations represent a pair of coupled oscillators, corresponding to the 
initial and the searched state, subject to a harmonic external field, with the 
initial conditions aj(0) = 1, an(0) = (n 7^ j). In such a case the last equation 
above leads to anit) = for all n s,j. The remaining two equations connect 
states s and j. Making the change of variables 
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aj{t) = Xj{t) exp [iVjj (costuo^ — 1) /^o] , 

as{t) = Xs{t) exp [iFss (cos^ot - 1) /uq] , (9) 



in those equations, we obtain 

dxj (t) 



dt 
dxs{t) 



—i fl Xs{t) exp [ie {coscuot — 1) /lvq] , 
—i Xj{t) exp [—is {cosujot — 1) /cuq] ■ (10) 



dt 

where e — Fss — Fjj. Uncouphng these equations we arrive at 
^^L^p. + is sinujot ^^j^^^ + Xj{t) = 0, 

CLv (JjL 

- tssmuot ^ + xs{t) = 0. (11) 
dt"^ dt 

In the exceptional case when £ = the external monochromatic field does not 
affect the functioning of the search algorithm. But, as a simple inspection of 
the differential equations (11) shows, when Fss 7^ Fjj and the ratio a = a;o/0 
is irrational their solutions are not periodic. This means that in general the 
application of the monochromatic field shall affect the functioning of the search 
algorithm. To be more quantitative, the last two equations, which differ only in 
the sign of the second term, may be transformed into Hill's - type differential 
equations after a change of variables x = exp {ie J sinu;o^ dt)y, where x is 
either Xj[t) or Xs{t). The solution of the resulting Hill's equations [11] is a 
linear combination of the functions 

yi(i) = exp(//i)ri(i), (12) 
y2(i) = exp(-//i)r2(i), 

where Ti(t), T2{t) are periodic functions with period tt, and the Floquet expo- 

l/i(^) + ^(^ 



nent fi obeys the relation cosh (vr/x) = ^ |/i(7r) + ^(tt) . The initial conditions 

for y,{t) and y2{t) are y^{0) = 1, 2/2(0) = 0; ^{0) = 0, ^^(0) = 1. The stabil- 
ity of the solutions (12) is determined by the values of the Floquet exponent 
/J,, which depends on the external field parameters. In particular, for zero or 
imaginary values of /i, the solutions are stable. 

Prom these solutions one may obtain the probabihty of the searched state 
Pg — \as(t)f. Fig. 1 maps the value of this probability, at the optimal time r, 
as a function of the parameters a and e. From it wc may conclude that the 
search algorithm is useful also in the presence of external field, since for an 
extense region of the field parameters it is possible to reach the searched state 
with high probability. 
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Fig. 1. Probability of finding the searched state at the optimal time as function of 
the field parameters a, e. 

4 Repeated measurements 



The search algorithm we are studying makes a transition from the initial to 
essentially the sought state, other states being negligibly populated. Thus, any 
measurement process will leave the system, with high probability, in one of 
these two states. The probability distributions associated to the searched state 
|s), and the initial state initially evolve according to the map of eq.(2). 
If at time ti the state is measured, the probabilities that the wavefunction 
collapses into either \s) or \j) are given by 

Pj(ti) ~ cos2(fi ti) , (13) 

respectively. If after this first measurement the system is in state |j) the 
probability distributions of the states \s) and \j) after a second measurement, 
at a time A^i after ti, are the eqs.(13), (the eqs.(13) exchanging s and j), cal- 
culated at Ati. Between consecutive measurements, the system undergoes an 
unitary evolution. Therefore, for arbitrary intervals Atj between consecutive 
measurements, the probability distributions of the states \s) and \ j) at a time 
ij+i — ti-\- Aij satisfy, always within the approximation ujnm ^ the matrix 
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equations 



Ps{ti+l) 



Pi Qi 
Qi Pi 



(14) 



where Pi = cos^(f2Atj) and qi — 1 — Pi. This equation is not a Markovian 
process because the transition probabilities are time-interval dependent. 

The general solution of the previous matrix equation for any time sequence of 
measurements is 



'Psm 

.PjiO))' 



(15) 



where P,(0) = 0, Pj{0) = 1 and 



ar 



1 
2 
1 
2 



n - 1] 

i=0 
m 

n [2p. - 1] 

1=0 



(16) 



If we now consider that the measurement processes are performed at regular 
time intervals, tn — nAi, eq. (16) becomes 



«m = 2 + (cos(2fiAt))' 
(5m = -[I- (cos(2QAi))' 



(17) 



Now, we have an unitary evolution between consecutive measurements, and 
the global evolution is a Markovian process, as all the Atj are equal. 

Let us now consider apphcations of the above eqs. (15-17). In the case where 
the value of the optimal time r is only approximately known, we take a single 
measurement at a time equal to our estimate of r, which we designate as r*. 
In that case the probability to find the searched state after the measurement 
is simply given by 



P,(r*) = ^[l-cos(2Qr*)]^l-^'^'^ 



r 



(18) 



which shows that the probability of finding the searched state, Ps{r*), de- 
creases quadratically with the relative error in the estimate of the optimal 
time. 



In Fig. 2 we present the probabilities of the searched and the initial states as 
a function of the number of measurements m, in a case where At = r*and the 
relative error is (r* — t)/t = 0.2. The full and the dashed hues correspond 
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Fig. 2. Probabilities of the searched (P^) and the initial (Pj) states as a function 
of the number of measurements (m), for fixed At = t*. The dashed and the full 
lines correspond to the theoretical development, eq. (17) in the same two cases, 
respectively. 

to the previous treatment, the circles (Pj) and the stars (P^) are the results 
obtained from a direct solution of the Schrodinger equation for the quantum 
rotor Hq, and for an ensemble of 500 trajectories with N — 50. We notice that 
after the first measurement the value of P, is close to 0.9. One also notices 
that after many measurements this value tends to 0.5. 

We remark that eq.(17) implies that for a sufficiently large number of mea- 
surements, both Pg and Pj tend to 1/2, independently of the interval between 
measurements At and the initial conditions. Further inspection of Fig. 2 shows 
that the average value of Pg over different number of measurements m is also 
1/2. This means that the algorithm can be still useful, as it may yield the 
searched state with a probability 0.5, even in the complete absence of knowl- 
edge on the number of states on which the search is performed. However, to 
perform this task in a reasonably small number of measurements m, one has 
to make sure At is not too small, see eq.(17). 

In the case where the m measurements are performed in a total time t = 
where the unperturbed algorithm converges to the searched state. At — r/m, 
so the coefficients am, Pm become 




(19) 
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Fig. 3. Probabilities of the searched (Pg) and the initial {Pj) states as a function of 
the number of measurements (m), performed in a total time r. The dashed and full 
lines are the theoretical results, obtained from eq.(19), in both cases, respectively. 

We show in Fig. 3 the probabilities Pg (stars), Pj (circles) after m measure- 
ments for = 50 and 500 trajectories as in the previous figure. The calculation 
based on the direct solution of the Schrodinger equation is compared to the 
one obtained from eq.(19), as a function of m. Besides the good agreement 
between the two calculations we notice that Pg decreases with m, so that for 
m > 30, Ps < 0.1. This simply means that the more frequently the wave 
function collapses, the harder it becomes for the algorithm to significantly de- 
part from the initial state. Therefore, in this case the algorithm behaves as an 
example of the Quantum Zeno effect, where the a high frequency of measure- 
ments hinders the departure of the system from its initial state [13,14]. This 
also explains our previous assertion, that for the algorithm to be useful the 
time At between measurements must not be too small. 



5 Conclusions 

We have extended the study of the search algorithm presented in [7]. There it 
was shown that the algorithm was robust when the energy of the searched state 
had some imprecision. In this work, the resonant algorithm was treated as an 
open system, and subject to two types of external interactions: an oscillating 
external field and measurement processes. 

It was shown that, although the algorithm is in general affected by a peri- 
odic external field, for extensive zones of the field parameter values it works 
with good efficiency. In the case of measurements, the probability distribution 
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for the searched and initial states were obtained analytically. When a set of 
periodic measurements are performed, the probability distribution satisfies a 
Master equation. In this way we have shown that the global behavior of the 
algorithm becomes a Markovian process. Therefore, the algorithm can still 
be useful in the case of repeated measurements, as long as they are not too 
frequent. 

We believe that these results may be directly extended to the Grover algo- 
rithm. As it may be seen from the description of this algorithm in ref. [5] , and 
as already pointed out by Grover [15], this algorithm is essentially a resonance 
between the searched and the average state. One should also be aware that, 
while in the case of the Grover algorithm care is needed to make sure that 
after each measurement the state of the system is reinitialized into the average 
state, this procedure is not required in the algorithm proposed in this work. 
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